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ABSTRACT 

We present multi- frequency radio continuum as well as H I observations of the compos- 
ite galaxy NGC 6764, which has a young, circumnuclear starburst and also harbours 
an active galactic nucleus (AGN). These observations have been made at a number 
of frequencies ranging from ~600 MHz to 15 GHz using both the Giant Metrewave 
Radio Telescope (GMRT) and the Very Large Array (VLA). They reveal the structure 
of the bipolar bubbles of non-thermal, radio emission which are along the minor axis 
of the galaxy and extend up to ~1.1 and 1.5 kpc on the northern and southern sides 
respectively. Features in the radio bubbles appear to overlap with filaments of Ha 
emission. The high-resolution observations reveal a compact source, likely to be asso- 
ciated with the nucleus of the galaxy, and a possible radio jet towards the south-west. 
We have compiled a representative sample of galaxies with bubbles of non-thermal 
radio emission and find that these are found in galaxies with an AGN. 

The Hi observations with the GMRT show two peaks of emission on both ends 
of the stellar-bar and depletion of Hi in the central region of the galaxy. We also 
detect H I in absorption against the central radio peak at the systemic velocity of the 
galaxy. The H i-absorption spectrum also suggests a possible weak absorption feature 
blue-shifted by ~120 km s _1 , which requires confirmation. A similar feature has also 
been reported from observations of CO in emission, suggesting that the circumnuclear 
starburst and nuclear activity affect the kinematics of the atomic and molecular gas 
components, in addition to the ionised gas seen in Ha and [Nil]. 

Key words: galaxies: individual: NGC6764 - galaxies: active - galaxies: starburst - 
galaxies: kinematics and dynamics - radio continuum: galaxies - radio lines: galaxies 



1 INTRODUCTION 

Galactic-scale outflows or galactic winds play an important 
role in heating and supplying kinetic energy to the inter- 
galactic medium (IGM) and enriching it with metals (e.g. 
Heckman, Armus & Miley 1990; Veilleux, Cecil & Bland- 
Hawthorn 2005 for a review). These outflows can be driven 
either by a circumnuclear starburst or an active galactic 
nucleus (AGN), although it is often likely to be due to a 
combination of both these processes. The conditions for an 
AGN, such as a deep potential well and a supply of gas, 
could also lead to a circumnuclear starburst and trigger a 
starburst-driven galactic wind. Observationally, about half 
of nearby Seyfert 2 galaxies are found to host a nuclear star- 
burst, with the proportion being larger for infra-red selected 
objects (e.g. Cid Fernandes et al. 2001; Veilleux 2001 and 
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references therein). Also in some cases, such as in NGC3079, 
there is strong evidence that starburst-driven winds co-exist 
with the AGN-driven outflows (Cecil et al. 2001; Irwin & 
Saikia 2003). Although it is usually not straightforward to 
disentangle the contributions to the galactic-scale outflows 
from the AGN and the starburst, the galactic winds from 
starbursts tend to be oriented along the minor axes of their 
parent galaxies, while the outflows from AGN show no pre- 
ferred orientation (e.g Ulvestad & Wilson 1984; Kinney et 
al. 2000; Gallimore et al. 2006 and references therein). 

In this paper, we present radio continuum and Hi obser- 
vations of the S-shaped spiral galaxy, NGC6764, which hosts 
an AGN and also exhibits evidence of a recent, intense star- 
burst. It has been classified as a Seyfert 2 galaxy (Rubin, 
Thonnard & Ford 1975) and also as a low-ionization nuclear 
emission line region or LINER (Osterbrock & Cohen 1982; 
Schinnerer, Eckart & Boiler 2000; Eckart et al. 1991, 1996). 
X-ray observations exhibit evidence of variability by a factor 
of ~2 over a time scale of 7 days, consistent with the exis- 
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Table 1. Basic data on NGC 6764. a 



RA (J2000) 6 


DEC (J2000) ft 


Type c 


a 


x b d 


v e 

V sys 


if 


log(L FJJ {)s 


D h 


(h m s) 


(°'") 




(' 


x ') 


(km s -1 ) 


(°) 


(L©) 


(Mpc) 


19 08 16.428 


50 55 59.47 


SB(s)bc 


2.3 


x 1.3 


2416 ±4 


44.5 


10.31 


34.0 



a Taken from the NASA Extragalactic Database (NED), unless stated otherwise. 

b The position of radio peak from our high-resolution, VLA A-array, 8460-MHz image. 

c Morphological type. 

d Optical major and minor axes. 

e Heliocentric systemic velocity. 

/ Inclination angle from Hi study by Wilcots, Turnbull & Brinks (2001). 

g Log of the far infra-red luminosity (Condon et al. 1996) revised for a distance of 34.0 Mpc. 

h Distance estimated using the galaxy recessional velocity and Ho=71 km s _1 Mpc -1 (Spergel et al. 2003) 

using the web-based cosmology calculator of Ned Wright. For this distance 1"=163 pc. 



tence of an AGN (Schinnerer et al. 2000). The basic proper- 
ties of the galaxy are summarised in Table 1. It is strongly 
barred and has a compact nuclear region, bright in both 
optical emission lines (Eckart et al. 1996) and radio contin- 
uum (Baum et al. 1993). There is a very prominent 466 nm 
Wolf-Rayet (W-R) emission feature from the nucleus (Os- 
terbrock & Cohen 1982; Schinnerer et al. 2000), suggestive 
of recent star formation (Armus, Heckman & Miley 1988; 
Conti 1991). There is a dense concentration of molecular 
gas in the central region (Eckart et al. 1991, 1996; Kohno et 
al. 2001) with Tkm > 20K in the circumnuclear region and 
<10K towards the spiral arms (Eckart et al. 1991). Interfer- 
ometric observations show molecular cloud complexes with 
an asymmetric velocity field which is blue-shifted relative to 
the systemic velocity. The distribution of the 2.12 /im H2 line 
flux is orthogonal to the stellar bar and also indicates pos- 
sible outflow from the circumnuclear region (Eckart et al. 
1996). The radio emission is extended roughly orthogonal 
to the major axis, forming bubbles of emission on opposite 
sides. There is also an extended X-ray component similar in 
extent to that of the radio emission (Schinnerer et al. 2000) . 

The new radio continuum observations presented here 
were made with the Very Large Array (VLA) and the Gi- 
ant Metrewave Radio Telescope (GMRT) at a number of 
frequencies ranging from ~600 MHz to 15 GHz. These ob- 
servations reveal the structure of the bubbles of non-thermal 
emission extending up to ~7 and 9 arcsec on the northern 
and southern sides respectively, as well as that of the circum- 
nuclear region, in far greater detail. We present spectral in- 
dex images of these features, compare NGC6764 with other 
sources with bubbles of non-thermal plasma and discuss the 
possible origin of these features. The Hi observations were 
made with the GMRT to examine the distribution of Hi 
clouds with higher resolution and also possible kinematic 
effects of the circumnuclear starburst on the Hi clouds via 
Hi-absorption spectra. We describe the observations in Sec- 
tion 2, and present the observational results from our radio 
continuum and Hi observations in Sections 3 and 4 respec- 
tively. These are discussed in Section 5 while the conclusions 
are summarised in Section 6. 



2 OBSERVATIONS AND DATA ANALYSIS 

The observing log for both the GMRT and VLA observa- 
tions is presented in Table 2, which is arranged as follows. 



Table 2. Observation log. 



Telescope 


Freq. 


Obs. 


t 


Phase 


Seal. 




MHz 


date 


min 


Calib. 


Jy 


GMRT 


598 


2004Jan04 


180 


2023+544 


0.80 


VLA-A 


1400 


2003Augl2 


27 


2023+544 


1.11 


GMRT 


1408,Hi 


2004Jan26 


240 


2023+544 


1.03 


VLA-A 


4860 


2003Augl2 


50 


2023+544 


1.08 


VLA-A 


8460 


2003Augl2 


60 


2023+544 


0.99 


VLA-A 


14965 


2003Augl2 


10 


2023+544 


0.87 


Archival 












VLA-B 1 


4985 


1985May25 


45 


1927+612 


0.55 


VLA-C 2 


8435 


1999July26 


16 


1944+548 


0.66 


VLA-B 3 


8460 


2002 


9 


1944+548 


0.65 


VLA-C 3 


14964 


2002 


16 


1824+568 


1.68 



Programme code: l : VH17; 2 : AL490; 3 : AC624 

Column 1 : Name of the telescope where we also list the con- 
figuration for the VLA observations. In addition to our own 
data, we have also analysed some of the relevant archival 
VLA data on this galaxy, which are listed separately in Ta- 
ble 2. Column 2: The frequency of the observations where H 1 
denotes spectral line observations. Columns 3 and 4: Dates 
of the observations and the time, t, spent on the source in 
minutes. For the archival data the dates have been taken 
from the image headers. However, for the observations in 
2002, the precise dates seem inconsistent with the array con- 
figurations; hence only the year is listed. The project codes 
for the archival data are indicated below the table. Columns 
5 and 6: The phase calibrator used and its flux density esti- 
mated from our observations. 

The observations with the GMRT, which is described on 
the website http://www.gmrt . ncra. tif r . res . in, as well as 
with the VLA were made in the standard fashion, with 
each source observation interspersed with observations of 
the phase calibrator. The primary flux density calibrator 
was 3C286 whose flux density was estimated on the Baars 
et al. (1977) scale, using the 1999.2 VLA values. The band- 
width of the continuum observations with the GMRT at 598 
MHz was 16 MHz, while for the L-band observations it was 
8 MHz. The bandwidth for all the VLA observations was 
50 MHz. The VLA data from both the IFs were combined 
to make the final images. The data analysis was done us- 
ing the Astronomical Image Processing System (AIPS) of 
the National Radio Astronomy Observatory. Since GMRT 
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NGC6764 OSS Fl-band 

^H^^^H^^^H^^^ NGC6764 NVSS 1400MHz 




Figure 1. Left panel: An optical R-band image of NGC 6764 from the Digitized Sky Survey (DSS). Contour levels are in arbitrary 
units. Middle panel: NVSS image of the same region with an angular resolution of 45 arcscc. Contour levels are 0.50x(— 4, —2.82, 2.82, 
4, 5.65, 8 . . . ) mjy/beam. All contour levels are in steps of \/2, unless otherwise specified. Right panel: GMRT 598-MHz image of the 
same region of the sky with an angular resolution of ~11 arcsec. Contour levels are 0.40x(— 4, —2.82, 2.82, 4, 5.65, 8 . . . ) mjy/beam. 
The cross denotes the position of the optical nucleus from Clements (1981) in all the images except Fig. 11, where the radio nucleus is 
marked. 
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Figure 2. The radio spectrum of the central region of NGC 6764. 
The dashed line represents the linear least-squares fit to the data 
points. 



data is acquired in the spectral-line mode with 128 spectral 
channels, gain and bandpass solutions were applied to each 
channel before combining them. 

The analysis of the H I observations was also done in 
the standard way. 3C286 was the primary flux density and 
bandpass calibrator. The total bandwidth for Hi observa- 
tions was 8 MHz and the spectral resolution was 62.5 kHz, 
which corresponds to 13.5 km s _1 in the centre of the band. 
We discarded a few antennas with more than 3% fluctuations 
in the bandpass gains during the observations. One channel 
in the beginning and eight channels towards the end were 
also not included in the analysis. The AIPS task UVLIN was 
used for continuum subtraction and the multi-channel data 
were then CLEANed using IMAGR with the robustness pa- 
rameter equal to zero. To estimate the overall properties of 
the galaxy, we made images by tapering the data to different 
uv distances. 

We have made self-calibrated images for all the differ- 
ent data sets. The GMRT images and the archival VLA B- 
and C-array 8460-MHz images showed significant improve- 



ment after self calibration. For these data sets two cycles of 
phase and one cycle of amplitude self calibration were ap- 
plied. For the remaining data sets, we have presented the 
un-selfcalibrated images. Some of the observed parameters 
of the GMRT and the VLA continuum images are presented 
in Table 3 which is arranged as follows. Columns 1 and 2 
are similar to that of Table 2, except that we also list the 
NRAO VLA Sky Survey (NVSS). Columns 3 to 5: The res- 
olution of the image with the major and minor axes being 
listed in arcsec and the position angle (PA) in degrees. Col- 
umn 6: The rms noise in units of mjy/beam. Columns 7 and 
8: The peak and total flux densities in units of mjy/beam 
and mjy respectively. These have been estimated by spec- 
ifying a polygon around the source. The total error in the 
flux density is approximately 5%. 

3 RADIO CONTINUUM EMISSION 
3.1 The radio bubbles 

In Fig. 1 we show the optical DSS R-band image of the 
galaxy (left panel), the NVSS image with an angular reso- 
lution of 45 arcsec (middle panel) and the GMRT 598-MHz 
contour map with an angular resolution of ~11 arcsec (right 
panel). The R-band image shows a strong stellar bar and 
peaks of emission which are likely to be Hn regions. The Ha 
continuum-subtracted image of the galaxy (Zurita, Rozas & 
Bcckman 2000) shows clearly the large number of Hn regions 
and also the filamentary-like structure extending beyond the 
circumnuclear region, defined here to be within ~0.5 kpc 
from the radio nucleus. At radio frequencies, the NVSS im- 
age appears somewhat extended with a deconvolved angular 
size of 15 x 14 arcsec 2 along a PA of 36°, but well- within the 
optical extent of the galaxy. The GMRT 598-MHz image 
shows the emission to be dominated by the central radio 
source, with possible weak emission from the bar. The spec- 
tral index of the central radio source (Fig. 2) which contains 
the bubbles of radio emission described below has been es- 
timated from our measurements to be — 0.74T0.02 (S (x i/ a ) 
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19 0817.5 17.0 16.5 16.0 15.5 

RIGHT ASCENSION (J2000) 



Figure 3. Left panel: VLA A-array image of the radio bubbles at 1400 MHz with an angular resolution of ~1.3 arcscc. Contour levels 
arc 0.052x(-4, -2.82, 2.82, 4, 5.65, 8 . . . ) mjy/beam in steps of %/2. Middle panel: VLA B-array image of the same region at 4985 MHz 
with an angular resolution of ~1.2 arcsec. Contour levels are 0.054x(— 4, —2.82, 2.82, 4, 5.65, 8 . . . ) mjy/beam. Right panel: Spectral 
index map of the region obtained by smoothing the 4985-MHz image to the same resolution as that of the 1400-MHz image. The contour 
level of —0.7 demarcates the regions of flatter and steeper spectral indices. 




Figure 4. Left panel: GMRT image of the radio bubbles in NGC 6764 at 1408 MHz with an angular resolution of ~2.7 arcscc. Contour 
levels are 0.15x(— 4, —2.82, 2.82, 4, 5.65, 8 . . . ) mjy/beam. Middle panel: VLA C-array image of this region at 8460 MHz with the same 
resolution as that of the GMRT image. Contour levels are 0.039x(— 4, —2.82, 2.82, 4, 5.65, 8 . . . ) mjy/beam. Right panel: Spectral index 
map of the region generated from these two images. The contour level of —0.7 demarcates the regions of flatter and steeper spectral 
indices. 



between between 598 and 8460 MHz. This shows the source 
to be clearly dominated by non-thermal emission. 

Published radio images of the central region of the 
galaxy (Wilson & Willis 1980; Ulvestad, Wilson & Sramek 
1981; Condon et al. 1982; Baum et al. 1993) have shown a 
'bright resolved peak at the location of the optical nucleus 
and diffuse emission on both sides of the nucleus' (cf. Baum 
et al. 1993). Our images of the extended emission with higher 



sensitivity using the VLA A-array at 1400 MHz and VLA 
B-array at 4985 MHz clearly show a bubble-like structure 
with a drop in brightness towards the centre of each bub- 
ble, located on opposite sides of the nucleus (Fig. 3: left and 
middle panels). These images have been made with angular 
resolutions of ~1.3 and 1.2 arcsec respectively. The extents 
of the bubbles are ~7 and 9 arcsec (1.1 and 1.5 kpc re- 
spectively) on the northern and southern sides respectively, 
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4860 MHz 
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NGC6764 
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Figure 5. Left panel: VLA A-array image of the circumnuclear region at 4860 MHz with an angular resolution of ~0.35 arcsec. Contour 
levels are 0.023X (—4, —2.82, 2.82, 4, 5.65, 8 . . . ) mjy/beam. Right panel: VLA A-array image of the circumnuclear region at 8460 MHz 
with an angular resolution of ~0.19 arcsec. Contour levels are 0.016x(— 4, —2.82, 2.82, 4, 5.65, 8 . . . ) mjy/bcam. 



while the flux densities at 1400 MHz are 42 mjy and 62 mjy 
respectively. In addition to the asymmetry in the extents of 
the bubbles along the minor axis of the galaxy, the radio 
emission in the circumnuclear region which is roughly or- 
thogonal to this axis is also asymmetric relative to the radio 
peak or the position of the optical nucleus, which is marked 
with a + sign. The emission on the eastern side extends for 
~5 arcsec (0.8 kpc) from the radio peak compared with 2.5 
arcsec (0.4 kpc) on the western side. At a somewhat lower 
resolution of 2.7 arcsec, our GMRT image at 1408 MHz and 
the VLA archival image at 8460 MHz with the same reso- 
lution (Fig. 4: left and middle panels) show the bubbles to 
have similar asymmetric morphologies. 

A spectral index image between 1400 and 4985 MHz 
(Fig. 3: right panel) made by smoothing the 4985-MHz im- 
age to the resolution of the 1400-MHz image and using only 
those pixels which are over 3 times the rms noise, shows that 
the mean spectral index is ~— 0.72. The western side appears 
to have a somewhat flatter spectral index, with a mean value 
of ~— 0.56, compared with ~— 0.77 for the rest of the source. 
The spectral index image between 1408 and 8460 MHz (Fig. 
4: right panel) shows a mean spectral index of ~— 0.85. This 
spectral index map also shows a similar trend of a flatter 
spectral index on the western side. The lateral or east-west 
asymmetry which is seen in the circumnuclear region is also 
apparent in parts of the bubbles, as seen in the VLA A-array 
image shown in Fig. 3. If this structural asymmetry, where 
the western side appears closer to the peak or the axes of the 
bubbles, and spectral asymmetry are related, it could be due 
to an asymmetric distribution of thermal gas. This gas could 
be probed by high-resolution X-ray observations. Although 
the ROSAT HRI image (Schinnerer et al. 2000) shows ex- 
tended X-ray emission in addition to the nuclear source, the 
resolution is inadequate. High-resolution X-ray observations 
with Chandra are required to study the detailed distribution 
of this gas. 

Adopting the spectral index of —0.74 as the mean value 
over the region of emission, and assuming the synchrotron 
emission to have lower and higher cut-offs at 10 7 and 10 10 



Hz, respectively, a proton-to-electron energy ratio of unity, a 
filling factor of unity and an oblate spheroidal distribution 
for the emitting region (17" x 9" x 9"), we estimate the 
minimum energy density and equipartition magnetic field to 
be 6.5 x 10~ 12 erg cm~ 3 and 8.4 /xG, respectively. The radia- 
tive lifetime of an electron radiating in this field at 1.4 GHz is 
2.5 xlO 7 yr. These values are similar to estimates for nearby 
galaxies (e.g. Condon 1992). The proton-to-electron ratio is 
not well determined for external galaxies. For a value of 50 
from studies of cosmic rays in our own Galaxy (e.g Webber 
1991), the equipartition magnetic field would increase by a 
factor of 2.5. 



3.2 The circumnuclear region 

The full-resolution VLA-A array images of the circumnu- 
clear region, which is taken to be within ~0.5 kpc from the 
radio nucleus, are presented in Fig. 5. These images at 4860 
and 8460 MHz with angular resolutions of ~0.35 and 0.19 
arcsec respectively reveal the structure to be somewhat com- 
plex. In the lower-resolution 4860-MHz image, the emission 
on the western side of the peak is more prominent and is 
triangular shaped, extending for about 3 arcsec (~ 0.5 kpc) 
towards the south and 2 arcsec (0.3 kpc) towards the west. 
The emission on the eastern side is more diffuse, faint and 
shorter in extent and contributes only 3.7 mjy of the total 
flux density of 17 mjy visible in the image. 

The higher-resolution 8460-MHz image reveals two ex- 
tensions from the radio peak, one along a PA of —90° and 
extending for ~1 arcsec (0.16 kpc) and the other which is 
more prominent and slightly longer extending for ~1.2 arc- 
sec (0.2 kpc) along a PA of -135°. Given Bridle & Perley's 
(1984) definition of a radio jet, the south-western feature 
could be classified as a possible jet which connects the radio 
nucleus with more diffuse emission which is oriented along 
a PA of ~180° and smoothly joins on to the ridge of emis- 
sion on the western side of the southern bubble (Figs. 5). 
However, this feature could also be a resolved structure at 
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Table 3. Observed parameters of radio continuum images. 



Telescope 


Freq. 


Beam size 




rms 


Spk 


Stot. 




MHz 


maj. 


min. 


PA 


mjy 


mjy 


mjy 






tt 


it 


O 


/b 


/b 




GMRT 


598 


1 O A 


O A 


33 


0.39 


136 


183 


NVSS 


1400 


45.0 


45.0 





0.50 


102 


113 


VLA-A 


1400 


1.33 


1.19 


7 


0.05 


14 


103 


GMRT 


1408 


2.89 


2.56 


-85 


0.14 


33 


109 


VLA-A 


4860 


0.37 


0.33 


17 


0.023 


1.6 


17 


VLA-A 


8460 


0.21 


0.18 


-16 


0.016 


0.8 


5.9 


VLA-A 


14965 


0.13 


0.11 


43 


0.221 


<0.7 


<0.7 


VLA-B 


4985 


1.31 


1.18 


-28 


0.054 


7.1 


43 


VLA-C 


8460 


2.89 


2.56 


-85 


0.039 


8.7 


25 


VLA-B 


8460 


1.53 


1.29 


-68 


0.069 


5.2 


20 


VLA-C 


14939 


2.32 


1.06 


-82 


0.166 


4.1 


14 



4 HI OBSERVATIONS 

Hi observations of NGC6764 have been reported earlier with 
an angular resolution of ~54 arcsec (8.8 kpc) using the VLA- 
D array (Wilcots, Turnbull & Brinks 2001). We observed this 
source with the GMRT with higher resolution to (i) iden- 
tify any kinematic effects of the AGN and starburst activity 
on the circumnuclear Hi gas via Hi 21-cm absorption lines 
towards the central continuum source, (ii) study the distri- 
bution of Hi gas relative to the Hn regions, and (iii) further 
investigate the nature of the observed depression in Hi col- 
umn density towards the centre of the galaxy. Wilcots et al. 
noted this depression but with the relatively coarse resolu- 
tion of their observations were unable to test the proposition 
that this might be due to absorption against the central con- 
tinuum source. The results of our observations, which have 
an rms noise of 0.6 mjy beam -1 per channel in the full- 
resolution images, are described below. 



'NG'C6764:'GMRT HI global profile' T - 
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Figure 6. GMRT Hl-cmission global profile of NGC 6764 mapped 
with a resolution of 37".76x 33". 36 along PA ~6° and smoothed 
to 27 km s _1 . 

the base of the bubble. The extensions to the east and west 
could also form the base of the bubbles seen in Fig. 3. 

3.2.1 The nuclear source 

The position of the radio peak determined from the VLA-A 
array 8460-MHz image is RA: 19 h 08 m 16. s 428, Dec 50° 55' 
59. "47 in J2000 co-ordinates (see Table 1). This is ~0.04 
arcsec from the peak in the VLA-A array 4860-MHz image 
(Fig. 5), but is ~0.55 arcsec from the position of the optical 
nucleus (RA 19 08 16.370; Dec: +50 55 59.58 in J2000 co- 
ordinates with an uncertainty of 0.22 arcsec) determined 
by Clements (1981). This difference between the radio and 
optical peaks may be due to higher extinction close to the 
nucleus of this galaxy. Similar differences have been seen in 
other galaxies as well, such as in NGC1482 (Hota & Saikia 
2005). Although the radio peak is likely to be associated 
with the AGN, it is undetected in our VLA-A array 14940- 
MHz image, the 3u upper limit being ~0.7 mjy. This is only 
marginally lower than the flux density of 0.8 mjy/beam at 
8460 MHz, suggesting that the spectrum could still be flat. 
A more sensitive VLA-A array 14940-MHz image is required 
to detect the possible nuclear source at higher frequencies 
and determine its spectrum. 



4.1 Hi emission 

In Fig. 6, we show the global profile obtained by tapering the 
data to 4 kA and making channel maps which have an rms 
noise of ~1 mjy/beam and an angular resolution of ~35 arc- 
sec (Fig. 7). In the channel maps, we confirm the extensions 
to the north (channel=55; 2515 km s" 1 ) and south (chan- 
nel=74; 2258 km s _1 ) noted by Wilcots et al., which are 
clearly seen in our data. The total Hi mass calculated from 
the global profile over the velocity range 2225 to 2562 km 
s _1 is 3.6 x 10 9 M Q . The total Hi mass detected by Wilcots 
et al. (2001) is (5.9±l)xl0 9 M Q (revised for a distance of 
34 Mpc), which implies that in our higher-resolution images 
we have resolved out some of the diffuse extended Hi-gas of 
the galaxy. Moment maps were generated by blanking each 
channel of this data cube at the 5<r level. The Hi column 
density image and the velocity field are shown in Fig. 8. 
The velocity field is similar to that published by Wilcots et 
al. (2001). The iso- velocity contours show evidence of non- 
circular motions due to the existence of a bar in NGC6764. 
We also see the velocity field being slightly distorted on the 
north-eastern and the south-western edges of the disk, as 
seen by Wilcots et al. (2001). The Moment-2 map shows the 
velocity width of the gas being slightly higher in the inner 
region than that in the outer region (Fig. 8). 

To study the Hi properties with higher resolution, we 
generated moment maps by tapering the data to 10 kA and 
blanking each channel of this data cube at the 5a level. The 
total Hi column density image with an angular resolution 
of ~15 arcsec shows the structure of the two well-resolved 
blobs of emission which are located on opposite sides of the 
circumnuclear region (Fig. 9: left panel). A superposition 
of the Hi column density image on the DSS R-band image 
of the optical galaxy (Fig. 9: middle panel) shows that the 
Hi line-emitting gas is located in the outer regions of the 
galaxy, and there appears to be no significant Hi emission 
towards the centre of the galaxy which harbours the young 
starburst and the AGN. This is unlikely to be due to absorp- 
tion against the central continuum source, which extends for 
only ^15 arcsec, similar to the resolution of these images, 
and contributes almost the entire continuum flux density 
seen in the NVSS image. The depression in the Hi column 
density in the central region extends over >50 arcsec. 
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Figure 7. NGC 6764 Hi-emission channel maps made with a resolution of 37". 76x33". 36 along PA ~6°. Channel 50 corresponds to a 
heliocentric velocity of 2602 km s , while channel 79 corresponds to a heliocentric velocity of 2211 km s . The systemic velocity of 
2416 km s~ 1 corresponds to channel 64. The velocity separation between adjacent channels is 13.48 km s . The crosses in all the images 
denote the position of the optical nucleus. Peak flux: 19.8 mjy/bcam. Contour levels are —4, —2.82, 2.82, 4, 5.65, 8 . . . mjy/beam. 



A comparison of the Hi column-density image with the 
Ha image of the galaxy (Zurita et al. 2000) shows that the 
most prominent region in Ha is the central starburst source 
which has no significant Hi gas. Although close to the cir- 
cumnuclear region this could be due to ionisation by the 
central starburst and AGN, the depletion of Hi at larger dis- 
tances could be due to an ISM phase transition from atomic 
to molecular hydrogen. As mentioned earlier, there is a con- 
centration of molecular gas extending over ~2.3 kpc along 
the bar of the galaxy (Eckart et al. 1991, 1996; Kohno et al. 
2001; Leon et al. 2003). 

The next most prominent regions of Ha emission are 
at the ends of the bar, with weaker emission along the bar 
and spiral arms. The peaks of Hi emission in the south- 



western blob appear displaced from the Ha peaks while in 
the northern lobe they appear to be roughly coincident at 
this resolution (Fig. 9; right panel). 



4.2 Hi absorption 

In Fig. 10, we present the Hi absorption profile towards the 
peak of the central continuum source (Fig. 4; left panel) with 
an angular resolution of ~2.7 arcsec and an rms noise of 0.6 
mjy beam -1 . No absorption is observed towards any other 
part of the radio source. The peak absorption occurs at a 
heliocentric velocity of 2426 km s _1 , which is consistent with 
the heliocentric systemic velocity 2416 km s _1 . In addition 
there is a weak blue-shifted component at a heliocentric ve- 
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Figure 8. Left panel: GMRT total intensity Hi contour as well as gray scale map of NGC 6764 made with a resolution of 37". 76x33". 36 
along PA~6°. The contour levels correspond to Hl-column densities of 5, 71 and 132 X 10 19 atoms cm -2 . Middle panel: The Hi- velocity 
field (Momcnt-1) generated from the same data set. The velocity contours from west to east are 2300, 2325, 2350, 2375, 2400, 2425, 2450, 
2475 and 2500 km s~ 1 . The systemic velocity is 2416 km s~ 1 . Right panel: Hl-velocity width map of NGC 6764 from the same data set. 
The contour levels are 13 and 26 km s _1 




Figure 9. Left panel: Hi contour map and gray-scale image of NGC 6764 made with a resolution of 17".37x 13". 77 along PA ~34°. The 
contour levels correspond to 1.21x(8, 11.31, 16, 22.62, 32..) in units of 10 20 atoms cm -2 . Middle panel: The Hi contour map with the 
same resolution superimposed on the DSS Red-band image. Right panel: The Hi contour map superimposed on the Ho image of Zurita 
ct al. (2000). 



locity of ~2300 km s _1 , which corresponds to a blue shift 
of ~120 km s _1 relative to the systemic velocity. Although 
this feature is very weak (~2cr) and requires confirmation 
from more sensitive observations, it is likely to be real since 
it is seen in both the Stokes RR and LL data sets and it 
extends over a few channels. It is interesting to note that 
mm-wavelength observations of CO have shown molecular 
gas in emission towards the circumnuclear region which is 
blue shifted by ~140 km s _1 relative to the systemic velocity 
(Leon et al. 2003). 

The peak optical depth, calculated using the back- 
ground peak continuum flux density of 33 mjy from the area 
of 9 arcsec 2 against the peak is 0.06. The total optical depth, 
calculated with data points over the full width at zero inten- 
sity (2184—2548 km s _1 ) is 0.5. Assuming a spin tempera- 
ture of 100 K this corresponds to a total Nhi ~ 1-23 x 10 21 
cm~ 2 . The total mass of the absorbing Hi clouds is 2.35 x 
10 6 M Q . The blue-shifted component spread over 2184-2346 
km s _1 has an optical depth of 0.13 which corresponds to 



an Nhi ~ 3.1 x 10 20 cm~ 2 . Although the blue shift could 
be due to non-circular motions due to the barred potential, 
such shifts may also be caused by gas clouds driven outwards 
by the circumnuclear starburst or AGN. In the latter sce- 
nario, the outflowing (~120 km s _1 ) Hi cloud, which has an 
estimated mass of 5.9xl0 5 M Q , assuming the cloud to have 
an area of 9 arcsec 2 , has a kinetic energy of 8.5 x 10 52 ergs. 



5 DISCUSSION 



5.1 Radio bubbles 

At radio continuum wavelengths, the most striking features 
are the bubbles of non-thermal plasma with a size of ~1.1 
and 1.5 kpc on the northern and southern sides respectively 
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Table 4. Comparision with other radio-bubble galaxies. 
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Figure 10. GMRT Hi-absorption spectrum taken from a 3"x3" 
region around the radio continuum peak of NGC 6764. The spatial 
resolution of the data cube is 2". 90x2". 55 along PA ~-87°; the 
spectrum has been smoothed to 27 km s _1 . 



of the circumnuclear region, and oriented roughly orthogo- 
nal to the stellar bar and the major axis of the galaxy (Fig. 
11). There also appears to be an extension on the eastern 
side extending beyond the base of the bubbles. In addition, 
smaller-scale structure with a size of ~0.5 kpc but asym- 
metric in extent relative to the radio nucleus is seen ori- 
ented roughly along the stellar bar. There is also a jet-like 
structure starting from the nucleus along an initial PA of 
— 135° which appears to connect the nucleus to more ex- 
tended emission towards the south-west. These features are 
shown in Fig. 11 where the smaller-scale structure is shown 
superimposed on the larger-scale bubbles. 

We first discuss the origin of these structures in this 
composite galaxy which harbours both a young starburst 
and an AGN. It is worth noting here that the activity in 
this galaxy does not appear to be triggered by interaction 
with any companion galaxy, since there is no evidence of any 
nearby galaxy within a distance of ~160 kpc (Wilcots et al. 
2001). However, the Hi observations do show evidence of a 



weak perturbation but this is unlikely to significantly affect 
the starburst and AGN activity seen in this galaxy. 

In order to understand the nature of these bubbles, we 
have compiled a representative sample of nearby galaxies 
with similar structures from the literature. A bubble is de- 
fined to be one whose width is larger than approximately 
half the size of its length. This working definition has been 
adopted to distinguish it from jets or jet-like features and 
normal lobes of radio emission formed by these jets whose 
axial ratios are usually higher than ~2 (e.g. Leahy, Muxlow 
& Stephens 1989). The sample of nearby galaxies with bub- 
bles is listed in Table 4, along with some of their properties. 
The table is arranged as follows. Column 1: name of the 
galaxy; column 2: galaxy classification from NED; column 
3: classification of the AGN from NED where Sy denotes a 
Seyfert galaxy, L a LINER galaxy, SB a starburst and Hn 
a galaxy with prominent Hn regions; column 4: distance to 
the galaxy in Mpc; columns 5 and 6: total extent of the 
bubbles in arcsec and kpc respectively. Columns 7 and 8: 
flux density and radio luminosity at 1.4 GHz of the bubbles 
and any core or jet emission. These have been estimated 
from the FIRST images, except for NGC2992 for which we 
have used the NED image (Ulvestad & Wilson 1989) and 
NGC6764 for which we have used the results presented in 
this paper. The flux density of the bubbles and any core 
or jet emission has been estimated by specifying a polygon 
around this region. The flux density value for the bubbles 
in Circinus is not available. Column 9: The FIR luminosity 
of the galaxies. Except for Circinus, for which we have esti- 
mated the value from Elmouttie et al. (1998b), these have 
been taken from Condon et al. (1990, 1998) and scaled to 
the distances listed in column 4. Column 10: References for 
the radio, X-ray and/or Ha images. 

It can be seen that the total extents of the bubbles range 
from ^0.4 to 15 kpc, with a median value of ~2.5 kpc, sim- 
ilar to that of NGC6764. While some bubbles lie within the 
scale-height of the interstellar medium (ISM) of the parent 
galaxy, others extend beyond it. In these radio bubbles, the 
ionised gas seen in optical emission lines such as Ha or [On] 
tend to be spatially related to the radio emission. The radio 
luminosity of the bubbles and any core or jet emission ranges 
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Figure 11. Contour map of the radio-continuum emission from 
the circumnuclear region of NGC6764 at 4860 MHz with the VLA 
A-array (Fig. 5, left panel) is superimposed on the gray scale 
image of the radio bubbles at 1400 MHz with the VLA A-array 
(Fig. 3, left panel). The bubbles, the eastern extension and a 
possible jet are labelled. The cross marks the position of the radio 
nucleus. The line at a PA of 73° shows the orientation of the 
stellar bar. 



from ' 

of ~1.4xl0 22 W Hz" 1 , while the infrared luminosity of the 
galaxies ranges from 9.47 to 10.98 in units of log(L F i R /L Q ) 
with a median value of ~ 10.3. Although the galaxies are 
luminous at infrared wavelengths, their luminosity is well 
below the threshold for defining these as ultra-luminous in- 
frared galaxies. 

In Table 4, it is striking that all the sources with non- 
thermal bubbles of radio emission have an AGN associated 
with it. It is worth comparing this with archetypal starburst 
galaxies such as M82 and NGC253, whose infrared luminosi- 
ties of 10.38 and 10.26 in units of log(L,FiR/Lo) respectively, 
are comparable with the sources in the sample of galaxies 
with bubbles. These two galaxies are dominated by the cir- 
cumnuclear starburst with no clear evidence of an AGN, al- 
though there have been attempts to identify one (Wills et al. 
1999; Mohan et al. 2002). High-resolution radio images of the 
nuclear regions of these two galaxies reveal a large number of 



compact radio components (Muxlow et al. 1994; McDonald 
et al. 2002 and references therein; Collison et al. 1994; Ul- 
vestad & Antonucci 1997 and references therein) which are a 
mixture of Hn regions and supernova remnants, while lower- 
resolution observations reveal a more extended halo of emis- 
sion (Seaquist & Odegard 1991; Reuter et al. 1992; Carilli et 
al. 1992). Although outflows along the minor axes are visible 
in both these two galaxies in X-rays and Ha (Strickland et 
al. 2004), there are no bubbles of non-thermal plasma similar 
to those seen in our sample of galaxies. The archetypal star- 
burst galaxy in the southern hemisphere, NGC1808, which 
is also dominated by a starburst with no unambiguous iden- 
tification of an AGN (Forbes et al. 1992; Jimenez-Bailon et 
al. 2005), has extended non-thermal emission with no clear 
bubbles of radio emission (e.g. Dahlem et al. 1990). The 
small-scale radio structure is dominated by a number of com- 
pact components (Saikia et al. 1990; Collison et al. 1994). 
Also, in a recent study of the superwind galaxy NGC1482, 
which has a remarkable hourglass-shaped optical emission- 
line outflow as well as bipolar soft X-ray bubbles of emission 
(Veilleux & Rupke 2002; Strickland et al. 2004), there is no 
non-thermal radio emission or bubbles associated with the 
superwind although radio emission is seen associated with 
the disk of the galaxy (Hota & Saikia 2005). Its infrared 
luminosity of log Lpir=10.66 is similar to the sources with 
bubbles listed in Table 4. Although the high-resolution ra- 
dio observations of the nuclear region of NGC1482 show a 
peak of emission with a steep radio spectrum, and more dif- 
fuse emission with secondary peaks, it is not clear if this 
feature could be associated with an AGN (Hota & Saikia 
2005). It is relevant to note here that, from optical spec- 
troscopic observations, Kewley et al. (2000) have classified 
it to be a starburst galaxy without an AGN. Considering 
these aspects, it is tempting to speculate that the formation 
of bubbles is closely linked to the existence of an AGN. It 
is interesting to note that Colbert et al. (1996) have high- 
lighted similar differences in the radio structures between 
Seyfert and starburst galaxies. 

Thus although the AGN possibly gives rise to the bub- 
bles of non-thermal radio emission, these would interact with 
the external environment in the host galaxy. A galactic wind 
generated by a starburst (cf. Veilleux et al. 2005) along with 
effects of buoyancy (Gull & Northover 1973; Stone, Wilson 
& Ward 1988) could lead to the bubbles of plasma being 
oriented along the minor axis of the galaxy. In this case, the 
radio jet close to the nucleus could be misaligned with the 
minor axis, as seen in NGC6764 as well as in the well-studied 
galaxy NGC3079 (cf. Irwin & Seaquist 1988; Kondratko et 
al. 2005). Using the formalism of Stone et al. (1988), we es- 
timate the time scale for the bubbles to reach the observed 
distances due to buoyancy to be ~35 Myr. This is likely to 
be an upper limit because of the presence of an AGN as 
well as galactic wind. The extended X-ray emission from a 
region similar to that of the radio bubbles (Schinnerer et al. 
2000) and the velocity structure in the Ha gas (Rubin et al. 
1975) suggests the presence of a galactic wind in this star- 
burst galaxy. Estimating an outflow velocity of ~80 km s" 1 
for the ionised gas observed by Rubin et al., the dynamical 
time scale for the Ha gas to reach distances of ~1.5 kpc is 
~20 Myr. For comparison, the circumnuclear star formation 
history suggests two starbursts with ages of 3— 5 Myr and 
between 15 and less than 50 Myr (Schinnerer et al. 2000), 
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Figure 13. The Ha and [Nil] velocity field measured along the 
minor-axis of the bubbles. Data plotted from the measurements 
of Rubin ct al. (1975). 




Figure 12. Top panel: Ha image of the circumnuclear region 
of NGC 6764 from the published image of Zurita ct al. (2000), 
showing the Ha-filaments. Bottom panel: Contour map of the 
radio bubble at 1400 MHz with the VLA A-array (Fig. 3, left 
panel) is superimposed on the same Ha image in gray scale. 



suggesting that both episodes of star formation could have 
affected the detailed structure of the bubbles. 



5.2 Kinematic effects on the ISM 

Interaction of radio jets and bubbles, as well as the galactic 
wind with the ambient medium could affect the kinematic 
properties of the ISM. However, in the case of NGC6764 
the radio features are likely to play a relatively smaller 
role. The total minimum energy in the radio features is 
only ~4.8xl0 54 ergs, which corresponds to a pressure of 
2.2xl0 -12 dynes cm' 2 , similar to values for the ISM in 
the Milky Way. The total energy in the radio features of 
NGC6764 are significantly smaller than typical values for 
the total energy in the outflows. Although the total energy 
in the various components in the outflow in NGC6764 is not 
known, the kinetic energy in the outflow of only the warm 
ionised and the neutral gas components in typical galaxies 
with galactic winds are ~10 53 to 10 56 and ~10 56 to 10 60 
ergs respectively (Veilleux et al. 2005). 

A continuum-subtracted Ha image of the central region 
of the galaxy (Zurita et al. 2000) shows a number of filamen- 
tary structures, reminescent of those seen in galaxies with 
galactic winds. We have determined the positional informa- 
tion of this Ha image by comparing the pixel positions of 
three stars which are visible in the Ha image with their po- 
sitions determined from the DSS image. A superposition of 
the radio continuum image at 1400 MHz with an angular 
resolution of 1.2 arcsec (0.2 kpc) on the Ha image (Fig. 12) 
shows some correlation of the radio emission with the Ha 
features, which are possibly driven outwards by a galactic 
wind. 

Spectroscopic observations of NGC6764 by Rubin, 
Thonnard & Ford (1975) showed that along the axis of the 
bubbles, the ionised gas appears red-shifted on the northern 
side and blue-shifted on the southern side by up to ~80 km 
s" 1 relative to the systemic velocity of the galaxy (Fig. 13). 
These motions are likely to be caused by a bi-polar outflow 
due to the circumnuclear starburst. 

In addition to the ionised gas component, the outflow 
due to the circumnuclear starburst could also affect the 
molecular and atomic gas components in the galaxy. For 
example, interaction of the bubbles with the molecular gas 
can give rise to CO gas compression and also affect their 
kinematic properties. Such effects have been seen in detailed 
studies of individual starburst galaxies. For example, in the 
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starburst galaxy NGC253, Sakamoto et al. (2006) find that 
the circumnuclear molecular disk harbouring the starburst is 
highly disturbed on small scales by individual young clusters 
and supernovae and globally by the well-known superwind. 
In the map of NGC6764 published by Eckart et al. (1991), 
the line ratio CO(2-1)/CO(1-0) shows higher values of ~0.9 
along the orientation of the bubbles but smaller values of 
^0.5 near the arms. Such higher values of the line ratio are 
expected when CO gas is compressed. Leon et al. (2003) find 
evidence of CO gas which is blue shifted by ~140 km s _1 and 
have a very high line ratio of ~2 for CO(2-1)/CO(1-0). As 
discussed earlier (Section 4.2) our Hi absorption spectrum 
shows a possible feature at a blue-shifted velocity of ~120 
km s~\ which is similar to that of the molecular gas ob- 
served in CO. A similar blue-shifted feature was also noted 
in the superwind galaxy NGC1482 by Hota & Saikia (2005), 
which was also interpreted to be due to the circumnuclear 
starburst. 

Hydrodynamical simulations of starburst-driven winds 
show that the cool gas could expand laterally in the disc 
of the galaxy, be carried vertically outwards by the tenuous 
superwind or be entrained in the interface between the hot, 
superwind fluid and the cool, dense ISM (e.g. Heckman et 
al. 1990; Heckman et al. 2000). The kinematic and physical 
properties of the different components discussed above are 
broadly consistent with the results of these simulations of 
winds. 



6 SUMMARY AND CONCLUDING REMARKS 

We have presented radio continuum and Hi observations of 
NGC6764, which has a young starburst as well as a Seyfert 
nucleus. The results are briefly summarised here. 

(i) The low-resolution GMRT image at 598 MHz with an 
angular resolution of ~11 arcsec (1.8 kpc) shows that almost 
all the radio emission from the galaxy comes from the central 
region of the source, with very little emission from the disk 
of the galaxy. The VLA A-array image at 1400 MHz with an 
angular resolution of ~1.3 arcsec accounts for more than 90 
per cent of the flux density visible in the NVSS image. This 
central source has a non-thermal spectrum with a spectral 
index of -0.74±0.02 between ~600 MHz and 8 GHz. 

(ii) With an intermediate resolution of 1.3 arcsec (0.2 
kpc) , the central source is resolved into two bubbles of non- 
thermal radio emission which are shell-like with a depression 
in the centre. The ridges of emission appear to overlap with 
Ha filaments seen in a continuum-subtracted Ha image of 
the galaxy. The bubbles are oriented along the minor axis 
of the galaxy with a total extent of ~2.6 kpc. The western 
edge of the bubbles tends to have a flatter spectral index 
than the rest of the bubbles. 

(iii) The image at 4860 MHz with a resolution of ~0.35 
arcsec (0.05 kpc) shows emission along the major axis of the 
galaxy as well as extensions which appear to connect to the 
ridges of emission in the bubbles. 

(iv) Our highest- resolution image (0.19 arcsec; 0.03 kpc) 
reveals a compact source, which is possibly associated with 
the nucleus of the galaxy, and is slightly offset (~0.55 arc- 
sec) from the position of the optical nucleus determined by 
Clements (1981). In addition there is a possible radio jet 



along a PA of —135° which appears to connect to more ex- 
tended emission towards the south. 

(v) A compilation of a representative sample of ten galax- 
ies with non-thermal bubbles shows that these are all associ- 
ated with an AGN. However, the detailed radio and optical 
structures, as well as the distribution of atomic and molec- 
ular gas, could be affected by a galactic wind due to a star- 
burst in addition to any AGN activity. Prominent starburst 
galaxies such as M82, NGC253, NGC1482 and NGC1808 
have similar infrared luminosity, but do not exhibit such 
bubbles of radio emission, suggesting that these are not usu- 
ally caused by the circumnuclear starburst. 

(vi) The H I observations, which are of higher resolution 
than those reported earlier by Wilcots et al. (2001), show 
that the two main peaks of emission are roughly coincident 
with the ends of the bar. There is a depletion of Hi towards 
the central region. Although this is likely to be due to the 
circumnuclear starburst and AGN close to the nucleus, at 
larger distances it could be due to phase transition in the 
ISM from atomic to molecular hydrogen. 

(vii) The peak absorption feature in the Hi-absorption 
profile has an optical depth of ~0.06 and a heliocentric veloc- 
ity of 2426 km s~\ consistent with the systemic velocity of 
the galaxy. The H I mass estimated from the absorption pro- 
file is ~2.4x 10 6 M . There is a suggestion of a weak absorp- 
tion feature at a blue-shifted velocity of 120 km s~\ which 
requires confirmation from more sensitive observations. A 
feature with a similar blue-shifted velocity has also been 
reported from observations of CO in emission, suggesting 
that the circumnuclear starburst and nuclear activity may 
have affected the kinematics of the atomic and molecular 
gas components of the ISM in the host galaxy, in addition 
to the ionised gas seen in Ha and [Nil]. 
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